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CI This paper aims at a comprehensive understanding on the novel elastic property of double- 

stranded DNA (dsDNA) discovered very recently through single-molecule manipulation techniques. 
A general elastic model for double-stranded biopolymers is proposed and a new structural parameter 
called the folding angle tp is introduced to characterize their deformations. The mechanical property 
of long dsDNA molecules is then studied based on this model, where the base-stacking interactions 
between DNA adjacent nucleotide basepairs, the steric effects of basepairs, and the electrostatic 
interactions along DNA backbones are taken into account. Quantitative results are obtained by 
using path integral method, and excellent agreement between theory and the observations reported 
by five major experimental groups are attained. The strong intensity of the base-stacking inter- 
actions ensures the structural stability of DNA, while the short-ranged nature of such interactions 
C/2 . makes externally-stimulated large structural fluctuations possible. The entropic elasticity, highly 

extensibility, and supercoiling property of DNA are all closely related to this account. The present 
work also suggests the possibility that negative torque can induce structural transitions in highly 
extended DNA from right-handed B-form to left-handed configurations similar with Z-form configu- 
ration. Some formulae concerned with the application of path integral method to polymeric systems 
are listed in the Appendix. 
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I. INTRODUCTION 

ly" - ) \ DNA molecule is the primary genetic material of most organisms. It is a double-helical biopolymer in which two 
■ chains of complementary nucleotides (the subunits whose sequence constitutes the genetic message) wind (usually 
right-handedly) around a common axis to form a double- helical structure jjj. Because of this unique structure, 
the clastic property of DNA molecule influences its biological functions greatly. There are mainly three kinds of 
deformations in DNA double-helix: stretching and bending of the molecule, twisting of one nucleotide chain relative 
to its counterpart. All these deformations have vital biological significance. During DNA replication, hydrogen bonds 
between the complementary DNA bases should be broken and the two nucleotide chains be separated. This strand- 
separation process requires cooperative unwinding of the double-helix [^J. In DNA recombination reaction, RecA 
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i-^j ■ proteins polymerize along DNA template and the DNA molecule is stretched to 1.5 times its relaxed contour length 
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It is suspected that thermal fluctuations of DNA central axis might be very important for RecA polymerization 



[5|. Another important example is the process of chromosome condensation during prophase of the cell cycle, where 
the long (circular) DNA chain wraps tightly onto histone proteins and is severely bent 0]. Further more, in living cells 
J> , DNA chain is usually closed, i.e., the two ends of the molecule is linked together by covalent bonds and the molecule 
becomes endless. With this chain-closing process, all those quantities characterizing the topological state of the chain 
are fixed and can only be changed externally by topoisomerases, which are capable of transiently cutting one or both 
DNA strands and making one strand pass through the other at the cutting point (in the case of type I topoisomerases 
||) or one segment of DNA pass through another (in the case of type II topoisomerases [0J§|). It is possible that 
this kind of enzyme-induced topology-changing processes are also closely related to the particular mechanic property 
of DNA molecule. For example, the frequency of collisions between two distant DNA segments in a circular DNA 
molecule is infl uence d by the different knot types and different linking numbers (for a definition of this quantity, see 



below and Sec. II C). A thorough investigation of the deformation and elasticity of DNA will enable us to gain better 
understanding on many important biological processes concerned with life and growth. 

Detailed study on DNA elasticity now becomes possible with the recent experimental developments, including, 
e.g., optical tweezer methods, atomic force microscopy, fluorescence microscopy. These techniques make it possible to 
manipulate directly single polymeric molecules and to record their elastic responses with high precision. Experiments 
done on double-stranded DNA (dsDNA) have revealed that this molecule has very novel elastic property P-[l6[ . When 
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a torsionally relaxed DNA is pulled with a force less than 10 picoNewton (pN), its elastic response can be quantitatively 
understood by regarding the chain as an inextensible thin string with certain bending rigidity (namely, the wormlike 
chain model |)|,[l7 18|). However, if the external force is increased up to 65 pN, DNA chain becomes highly extensible. 
At this force, the molecule transit to an over-stretched configuration termed S-DNA, which is 1.6 times longer than the 
same molecule in its standard B-form structure [0,0. Besides external forces, it is also possible to apply torsional 
constraints to DNA double-helix by external torques. The linking number of DNA, i.e., the total topological turns 
one DNA strand winds around the other, can be fixed at a value larger (less) than the molecule's relaxed value. In 
such cases we say the DNA molecule is positively (negatively) supercoiled. It is shown experimentally ]p| that when 
external force is less than a threshold value of about 0.3 pN, the extension of DNA molecule decreases with increasing 
twisted stress and the elastic response of positively supercoiled DNA is similar to that of negatively supercoiled DNA, 
indicating the DNA chain might be regarded as achiral. However, if the external force is increased to be larger than 
this threshold, negatively and positively supercoiled DNA molecules behave quite differently. Under the condition 
of fixed external force between 0.3 pN and 3 pN, while positive twist stress keeps shrinking the DNA polymer, the 
extension of negatively supercoiled DNA is insensitive to supercoiling degree [l^Ji^] . In higher force region, it is 
suggested by some authors that positivel y s upercoiled DNA may transit to a configuration called Pauling-like DNA 
(P-DNA) with exposed nucleotide bases Jig] , while negative torque may lead to strand-separation in DNA molecule 
(denaturation of DNA double- helix jjj} ) . A very recent systematic observation performed by Leger et al. jl6| , on the 
other hand, suggested another possibility that negative supercoiling may result in left-handed Z-form configuration 
in DNA. 

The above-mentioned complicated elastic property revealed by experiments may be directly related to the versatile 
roles played by DNA molecule in living organisms. Theoretically, to understand DNA elastic property is of current 
interest. Concerned with one or another aspect of DNA elasticity, models were proposed and valuable insights were 
obtained (see, for example, Refs. [l7 28 1), and now it is widely accepted that the competition between DNA bending 
and torsional deformations deserves to pay considerable attention in understanding the elastic property of DNA 
molecules. However, it is still a great challenge to understand systematically and quantitatively all aspects of DNA 
mechanical property based on the same unified framework. What is the intrinsic reason for DNA molecule's entropic 
elasticity, highly extensibility as well as its supercoiling property? Is it possible for negative torque to stabilize 
left-handed DNA configurations? These are just some examples of unsolved questions. 

In the present work, we have tried to obtain a comprehensive and quantitative understanding on DNA mechanical 
property. We have thought that the double-stranded nature of DNA structure should be extremely important to its 
elastic property, and therefore have constructed a general elastic model in which this characteristic is properly taken 
into account via the introduction of a new structural parameter, the folding angle ip. The elastic property of long 
dsDNA molecules was then studied based on this model, where the base-stacking interactions between DNA adjacent 
nucleotide basepairs, their steric effects, and the electrostatic interactions along DNA backbones were all considered. 
Quantitative results were obtained by using path integral method, and excellent agreement between theory and the 
experimental observations of several groups were attained. It was revealed that, on one hand, the strong intensity of 
the base-stacking interactions ensures the structural stability of DNA molecule; while on the other hand, the short- 
ranged nature of such interactions makes externally-stimulated large structural fluctuations possible. The entropic 
elasticity, highly extensibility, and supercoiling property of dsDNA molecule are all closely related to this fact. The 
present work also revealed the possibility that negative torque can induce structural transitions in highly extended 
DNA from right-handed B-form configuration to left-handed Z- form-like configurations. Some discussions on this 
respect were performed and we suggested that a possible direct way to check the validity of this opinion is to measure 
the values of the critical torques under which such transitions are anticipated to take place by the present calculations. 

This paper is organized as follows: In Sec. [n] we introduce the elastic model for dsDNA biopolymers. At the 
action of an external force, the elastic response of dsDNA molecules is investigated in Sec. [II and compared with 
experimental observations of Smith et al. [p||l2"[ and Cluzel et al. [jllj. Section |y| focuses on supercoiled dsDNA 
molecules, where the relationship between extension and supercoiling degree is obtained numerically and compared 
with the experiment of Strick et al. E^] . From the calculated folding angle distribution, we infer that negative torque 
can cause structural transitions in dsDNA molecules from right-handed double-helix to left-handed ones. Section 
is reserved for conclusion. Two appendices are also presented: in Appendix [X] we review some basic ideas on 
the application of path integral method in polymer physics; and in Appendix [b| we list the matrix elements of the 
operators in Eq. ( |l5|) and Eq. (23). Some parts of this work have been briefly reported in a previous letter pg| . 
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II. ELASTIC MODEL OF DOUBLE-STRANDED DNA MOLECULE 



As already stressed, DNA molecule is a double-stranded biopolymer. Its two complementary sugar-phosphate chains 
twist around each other to form a right-handed double-helix. Each chain is a linear polynucleotide consisting of the 
following four bases: two purines (A, G) and two pyrimidines (C, T) pLpj. The two chains are joined together by 
hydrogen bonds between pairs of nucleotides A-T and G-C. Hereafter, we refer to the two sugar-phosphate chains as 
the backbones and the hydrogen-bonded pairs of nucleotides as the basepairs. In this section we discuss the energetics 
of such an elastic system (see Fig. [l]). First of all, the bending energy of the backbones of such double-stranded 
polymers will be considered; then we will discuss the interactions between DNA basepairs and energy terms related 
with external fields. 



A. Bending and folding deformations 



The backbones can be regarded as two inextensible wormlike chains characterized with a very small bending rigidity 
k = ksT£ p , where ks is Boltzmann's constant and T the environmental temperature, and l p ~ 1.5 nm is the bending 
persistence length of single-stranded DNA (ssDNA) chains jjJI . The bending energy of each backbone is thus expressed 
as k Jq '(dti/ds) 2 ds, where tj(s) (i = 1,2) is the unit tangent vector at arclength s along the i-th backbone pfj| , ^l| , 
and L is the total contour length of each backbone. The position vectors of the two backbones are expressed as 
r i (s) = f s U(s')ds'. 

Since there are many relatively rigid basepairs between the two backbones in many cases the lateral distance 
between the backbones can be regarded to be constant and equal to 2i?||] In this subsection we focus on the bending 
energy of the backbones, therefore, for the moment we regard each basepair as a rigid rod of length 2R linking 
between the two backbones and pointing along direction denoted by a unit vector b from iq to r 2 (Fig. [I]). Then, 
r 2 (s) — i"i(s) = 2i?b(s). In B-form DNA the basepair plane is perpendicular to DNA axis, therefore, in our model 
relative sliding of the two backbones is not considered and the basepair rod is thought to be perpendicular to both 
backbones pl| , with b(s) • ti(s) = b(s) • ta(s) = 0. The central axis of the double-stranded polymer can be defined 
as r(s) = iq(s) + Rh(s) [= r 2 (s) — Rh(s) = (ri(s) + r 2 (s))/2], and its tangent vector is denoted by t. In consistence 
with actual DNA structures, the central axial tangent t is also perpendicular to b, i.e., b(s) • t(s) = 0. (Notice that, 
however, t(s) ^ dr/ds; in this paper, s always refers to the arclength of the backbones.) 

Since all the tangent vectors ti, t 2 , and t lie on the same plane perpendicular to b, we can write that 



ti(s) = t(s) cos<p(s) + n(s) sinews), 
t 2 (s) = t(s) cos^(s) — n(s) sinews), 



(1) 



where n is also a unit vector and n = b x t, and ip is defined as half the rotational angle from t 2 to ti, with b 
being the rotational axis (Fig. |l|). We call ip the folding angle, and it can vary in the range (— 7r/2, +7r/2), with 
Lp > corresponding to right-handed rotations and hence right-handed double-helical configurations and tp < to 
left-handed ones. With the help of Eq. (@), we know that 

db t 2 -ti sin ip 

and 

dr 1 . 
— = -(ti + t 2 ) = t cost/?. (3) 
ds 2 

Equation indicates that cos (p measures the extent to which the backbones are "folded" with respect to the central 
axis. Based on Eqs. (]l]-0), the total bending energy of the two backbones can be expressed in the following form: 



In our present work, we have not taken into account the possible deformations of the nucleotide basepairs. In many cases 
this may be a reasonable assumption. However, under some extreme conditions such kind of deformations may turn to be 
important. For example, when DNA double-helix is stretched and at the same time is applied with a large positive torque, the 
nucleotide basepairs may collapse (see Ref. [H for a description). 



3 




The bending energy is thus decomposed into the bending energy of the central axis (the first term of Eq. (Q)) plus the 
folding energy of the backbones (the second and third terms of Eq. (§])). The physical meanings of these two energy 
contributions are very clear, and Eq. (^) is very helpful for our following calculations. In Eq. the bending energy 
of the central axis is very similar with that of a wormlike chain |l8| , both of which are related to the square of the 
changing rate of the axial tangent vectors. But there are two important differences: (a) in the derivative dt/ds of Eq. 
(||), the arclength parameter s is measured along the backbone, not along the central axis; (b) in the wormlike chain 
model the central axis is inextensible, while here the central axis is extensible. 




Backbones 

Central axis 

Basepairs 




FIG. 1. Schematic representation of a double-stranded DNA model used in this paper. The right part demonstrates the 
definition of 95 on the local t — n plane, where t, ti and t2 are, respectively, the tangential vectors of the central axis and the 
two backbones; is the folding angle; and the unit vector b = t x n is perpendicular to the t — n plane. 

In deriving Eq. the basepairs are models just as thin rigid rods of fixed length, i.e., the DNA molecule 

is viewed as a ladder-like structure (see Fig. |l|). Actually, however, basepairs form disc- like structures and have 
finite volume. The steric effects caused by the finite volume of basepairs is anticipated to hinder considerably the 
bending deformation of the central axis, and hence will increases its bending rigidity greatly Q. Furthermore, dsDNA 
molecule is a strong polyelectrolyte, with negatively-charged groups distributed regularly along the chain's surface. 
The electrostatic repulsion force between these negatively-charged groups will also considerably increase the bending 
rigidity of the dsDNA chain pL To quantitatively take into account the above mentioned two kinds of effects is 
very difficult. Here we treat this problem phenomenologically by simply replacing the bending rigidity k in the first 
term of Eq. (0) with a quantity n*. It is required that n* > k, and the precise value of k* will then be determined 



self-consistently by the best fitting with experimental data as shown in Sec. [II. 



B. Base-stacking interactions between basepairs 



In Sec. II A, we have discussed in detail the bending energy of dsDNA polymers, which is caused by bending 
of the backbones as well as steric effects and electrostatic interactions. In dsDNA molecule there is another kind 
of important interactions, namely the base-stacking interaction between adjacent nucleotide basepairs The 
base-stacking interactions originate from the weak van der Waals attraction between the polar groups in adjacent 
nucleotide basepairs. Such interactions are short-ranged and their total effect is usually described by a potential 
energy of the Lennard- Jones form (6-12 potential Q). Base-stacking interactions play significant role in stabilization 
of DNA double-helix. The main reason why DNA can but RNA can not form long double-helix is as follows |32] : 
Because of the steric interference caused by the hydroxyl group attached to the 2' carbon of RNA riboses, the stacking 
interaction between adjacent RNA nucleotide basepairs is very weak and can not stabilize the formed double-helical 
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structure; while in the DNA ribose, it is a hydrogen atom attached to its 2' carbon and serious steric interference is 
avoided (fortunately!). 

In a continuum theory of elasticity, the summed total base-stacking potential energy is converted into the form of 
the following integration: 



N-i . L 
Elj = y~] Ui, i+ i = / p(ip)ds, 
i=i Ja 



(5) 



where £/i,j+i is the base-stacking potential between the i-th and the (i + l)-th basepair, N is the total number of 
basepairs, and the base-stacking energy density p is expressed as 

~ I >o - 2cosVo] (for^<0). (bj 

In Eq. the parameter ro is the backbone arclength between adjacent bases (ro = L/N); (po is a parameter related 
to the equilibrium distance between a DNA dimer (rocos^o ~ 3.4 A); and e is the base-stacking intensity which is 
generally base-sequence specific . In this paper we focus on macroscopic properties of long DNA chains composed 
of relatively random sequences, therefore we just consider e in the average sense and take it as a constant, with e 
~ 14.0 ksT as averaged over quantum-mechanically calculated results on all the different DNA dimers Q. 

The asymmetric base-stacking potential Eq. (^J) ensures a relaxed DNA to take on a right-handed double-helix 
configuration (i.e., the B-form) with its folding angle <p ~ (fo. To deviate the local configuration of DNA considerably 
from its B-form generally requires a free energy of the order of e per basepair. Thus, DNA molecule will be very 
stable under normal physiological conditions and thermal energy can only make it fluctuate very slightly around its 
equilibrium configuration, since e 3> kT . Nevertheless, although the stacking intensity e in dsDNA is very strong 
compared with thermal energy, the base-stacking interaction by its nature is short-ranged and hence sensitive to the 
distance between the adjacent basepairs. If dsDNA chain is stretched by large external forces, which cause the average 
inter-basepair distance to exceed some threshold value determined intrinsically by the molecule, the restoring force 
provided by the base-stacking interactions will no longer be able to offset the external forces. Consequently, it will 
be possible that the B-form configuration of dsDNA will collapse and the chain will turn to be highly extensible. 
Thus, on one hand, the strong base-stacking interaction ensures the standard B-form configuration to be very stable 
upon thermal fluctuations and small external forces (this is required for the biological functions of DNA molecule 
to be properly fulfilled 0); but on the other hand, its short-rangedness gives it considerable latitude to change its 
configuration to adapt to possible severe environments (otherwise, the chain may be pulled break by external forces, 
for example, during DNA segregation ). This property of DNA base-stacking interactions is very important to 
dsDNA molecule. As we will see in Sees. [II and |y[ the mechanical property of DNA chain is indeed closely related 
to the above-mentioned insight. 



C. External forces and torques 

In the previous two subsections, we have described the intrinsic energy of DNA double-helix. Experimentally, to 
probe the elastic response of linear DNA molecule, the polymer chain is often pulled by external force fields and/or 
untwisted or overtwisted by external torques. To study the mechanic response of dsDNA molecule, we consider in 
this subsection the energy terms related to external forces and torques in our theoretical framework. 

For the external force fields, here we constrained ourselves to the simplest situation where one terminal of DNA 
molecule is fixed and the other terminal is pulled with a force F = Jzq along direction of unit vector zq ||. (In fact, 
hydrodynamic fields or electric fields are also frequently used to stretch semiflexible polymers Q, but we will not 
discuss such cases in this paper.) The end-to-end vector of a DNA chain is expressed as t(s) cosip(s)ds according 
to Eq. (|^) . Then the total "potential" energy of the chain in the external force field is 

Ef — — t cos ipds ■ F = — / ft ■ zq cos tpds. (7) 
Jo Jo 

In the experimental setup, external torques can be applied on linear dsDNA molecule by the following procedure: 
first, DNA ligases are used to ligate all the possible single-stranded nicks; then the two strands of dsDNA molecule 
at one end are fixed onto a template, while the two strands at the other end is attached tightly to a magnetic bead; 
afterwards, torques are introduced into DNA double-helix by rotating the magnet bead with an external magnetic field 
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PlP^f , Torque energy is then related to the topological turns caused by the external torque on DNA double- helix. 
The total number of topological turns one DNA strand winds around the other, which is usually termed the total 
linking number Lk |34-36|, is expressed as the sum of the twisting number, Tw(t\{t2), r) of backbone ri (or i^) 



around the central axis r and the writhing number Wr(r) of the central axis; i.e., Lk = Tw + Wr. According to Refs. 
13-f37| and Eq. (§), w <3 obtain that 



Zir 



L , db , 1 

t x b • —as = — 



ds 



2tt 



sin tp 
R 



ds. 



(8) 



The writhing number of the central axis is generally much more difficult to calculate. It is expressed as the following 
Gauss integral over the central axis Q]: 



Wr(r) 



1 

47T 



dr x dr' ■ (r — r' 

It- _ v'|3 



(9) 



In the case of linear chains, provided that some fixed direction (for example the direction of the external force, Zq) 
can be specified and that the tangent vector t never points to — zq (i.e., t - zq 7^ —1), it was proved by Fuller that the 
writhing number Eq. IM) can be calculated alternatively according to the following formula p8[ : 



r \ 1 /* z o x * • dt /ds , 

Wr(r) = — / ■ — ds. 

w 2tt J 1 + z • t 



(10) 



The above equation can be further simplified for highly extended linear DNA chains whose tangent t fluctuates only 
slightly around z . In this case, Eq. (nCl) leads to the approximate expression that 



Wr{r) 



4tt .,„ 



dt x 
1 ds 



dty 

ds 



ds. 



(11) 



where t x and t y are, respectively, the two components of t with respect to two arbitrarily chosen orthonormal directions 
(xo and yo) on the plane perpendicular to z . 

The energy caused by the external torque of magnitude T is then equal to 



Ei 



-2irTLk = -2itY{Tw + Wr). 



(12) 



To conclude this section, the total energy of a dsDNA molecule under the action of an external force and an external 
torque is expressed as 



E — Ef, + El j + Ef + E t 



ds ds 



— sin 4 tp + p(ip) — ft ■ zq cos tp — — sin tp — T 



zq x t ■ dt/ds 



R 2 



R 



1 + z • t 



ds 



— sin tp + p(tp) - ft ■ z cos^ - —smtp+ -t y - 77*3 j 



R 2 



R 



ds 



ds 



ds. 



(13) 
(14) 



Notice that Eq. (^) can be applied only in the case of highly extended DNA. In the following two sections, we will 
study the mechanical property of single dsDNA molecules based on the model energy Eq. (|l3|) and Eq. (14). The 
theoretical results will be compared with experimental observations and discussed. 



III. EXTENSIBILITY AND ENTROPIC ELASTICITY OF DNA 



In this section we investigate the elastic responses of single DNA molecules under the actions of external forces 
based on the model introduced in Sec. ||. There is no external torque acted, thus T = in Eq. (13). The particular 
form of the energy function Eq. ([l3]) of the present model makes it convenient for us to study its statistical property 
by path integral method. In Appendix ^ a detailed description on the application of path integral method to polymer 
physics is given p9[ |. Our calculations in this and the next sections are based on this method. 

For a polymer whose energy is expressed in the form of Eq. ( |i"3| ) with T — 0, according to the technique outlined in 
Appendix |A| (see Eqs. (Al) and (|A7|) ) , the Green equation governing the evolution of the "wave function" ^(t,tp;s) 
of the system is obtained to be of the following form: 



G 



ds 



t) 2 



d 2 



Al* p dt 2 M p dcp 2 k B T 



fcosip p(ip) 

t ■ Zq - 



k B T R 2 



V ■ 4 

— sin tf 



(15) 



where I* — k* /k B T and I = n/k B T. T he sp ectrum of the above Green equation is discrete and, for a long dsDNA 
molecule according to Eqs. (A10) and (A13), its average extension can be obtained either by differentiation of the 
ground-state eigenvalue, go, of Eq. jl5| ) with respect to /: 



(Z) = [ (t • z cos(p)ds = -Lk B T^j, 
Jo °7 

or by a direct integration with the normalized ground-state eigenfunction, $o(t, <f), of Eq. (p.5| 

(Z) = L J <f> | 2 t • z cos tpdtdip. 



(16) 



(17) 



Both go an < l > o(t, <p) can be obtained numerically through standard diagonalization methods and identical results are 
obtained by Eqs. ( plf ) and (17). Here we just briefly outline the main procedures in converting Eq. ( [l5| ) into the form 
of a matrix. 




FIG. 2. Force-extension relation of torsionally relaxed DNA molecule. Experimental data is from Fig. 2A of (symbols). 
Theoretical curve is obtained by the following considerations: (a) £ p = 1.5 nm and e = 14.0fc.eT; (b) £* = 53.0/2(cos <p)f=o 
nm, ro = 0.34/ (cos <p)f=o nm and R = (0.34 x 10.5/27r)(tany)/=o nm; (c) adjust the value of ipo to fit the data. For each ipo, 
the value of (cos <p) /=o is obtained self-consistently. The present curve is drawn with ipo = 62.0° (in close consistence with the 
structural property of DNA), and (cosy>)/=o is determined to be 0.573840. DNA extension is scaled with its B-form contour 
length L{cosip) f—o. 

Firstly, for our convenience we perform the following transformation: 

V = V-|, (18) 

hence the new argument ip can change in the range from to n. Then, we choose the combination of Y/ m (t) and 
f n {0) as the base functions of the Green equation Eq. (|l5|): 

*(t,£;«) =J2 C imn(s)Yi m (t)f n (0). (19) 
hnn 
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In the above expression, Y; m (t) = Yi m {6, 4>) (I = 0, 1, 2, ■ • • ; m = 0, ±1, • ■ • , ±Z) are the spherical harmonics pQ] , where 
6* and are the two directional angles oft, i.e., t = (sin cos sin 9 sin 0, cos 9); and 

/2 ri7r 
-sin(— <p) (n = l,2,..-) (20) 
a a 

are the eigenf unctions of one-dimensional infinitely deep square potential well of width a. f\ 




FIG. 3. Folding angle distribution for torsionally relaxed DNA molecules under external forces. 

With the wave function 'F being expanded using the above mentioned base functions, the operator acting on 'F (i.e., 
the expression in the square brackets of Eq. (^)) can also be written into matrix form under these base functions. 
This matrix, whose elements are listed in Appendix is then diagolized numerically to obtain its ground-state 
eigenvalue and eigenfunction. To simplify the calculation, we further notice that in the present case of Eq. (|J), the 
ground-state is independent to 4>, i.e., m can be set to m = in Eq. (|ig|). 

The resulting force vs extension relation obtained from Eq. (|l6|) or Eq. (|l7]) is shown in Fig. ^ in the whole 
relevant force range and compared with the experimental observation of Cluzel et al. [jll 12 1. The theoretical curve 
in this figure is obtained with just one adjustable parameter (see caption of Fig. |^); the agreement with experiment 
is excellent. Figure || demonstrates that the highly extensibility of DNA molecule under large external forces can be 
quantitatively explained by the present model. 

To further understand the force-induced extensibility of DNA, in Fig. || the folding angle distribution of dsDNA 



molecule is shown, with the external force kept at different values. Here, according to Eq. (A13) the folding angle 
distribution P(<p) is calculated by the following formula: 

P(<p) = [ |$ (t,^)| 2 dt. (21) 



2 In the actual calculations, the right boundary of the square well is chosen to be slightly less than it to avoid flush-off of 
computer memory caused by the divergence of the base-stacking potential Eq. Q at <p — tv. We set a — 0.95-7T in this paper. 
However, we have checked that the results are almost identical for other values of a, provided that a > 165°. 
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Figures || and 0, taking together, demonstrate that the elastic behaviors of dsDNA molecule are radically different 
under the condition of low and large applied forces. In the following, we will discuss them separately. 
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FIG. 4. Low-force elastic behavior of DNA. Here experimental data is from Fig. 5B of ]9j, the dotted curve is obtained for a 
wormlike chain with bending persistence length 53.0 nm and the parameters for the slid curve are the same as those in Fig. ^. 

The low- force region When external force is low (< 10 pN), the folding angle is distributed narrowly around 
the angle of ip ~ +57°, and there is no probability for the folding angle to take on values less than 0° (Fig. ||), 
indicating that DNA chain is completely in the right-handed B-form configuration with small axial fluctuations. This 
should be attributed to the strong base-stacking intensity, as pointed out in Sec. [IB. Consequently, the elasticity 



of DNA is solely caused by thermal fluctuations in the axial tangent t (Fig. 0), and DNA molecule can be regarded 
as an inextensible chain. This is the physical rc^son_ why, in this force region, the clastic behavior of DNA can be 
well described by the wormlike chain model |T~ 
chain model and the present model give identica 



18 41 . Indeed, as shown in Fig. [|, at forces < 10 pN, the wormlike 
results. Thus, we can conclude with confidence that, when external 



fields are not strong, the wormlike chain model is a good approximation of the present model to describe the elastic 
property of dsDNA molecules; and the bending persistence length of the molecule is 2£* (cos <p), as indicated by Eq. 
© and Eq. |). 

The large-force region With the continuous increase of external pulling forces, the axial fluctuations becomes 
more and more significant. For example, at forces ~ 50 pN, although the folding angle distribution is still peaked at 
ip ~ 57°, there is also considerable probability for the folding angle to be distributed in the region ip ~ 0° (Fig. ||). 
Therefore, at this force region, DNA polymer can no longer be regarded as inextensible. At / ~ 65 pN, another peak 
in the folding angle distribution begins to emerge at ip ~ 0°, marking the onset of cooperative transition from B-form 
DNA to overstretched S-form DNA 



I! 



This is closely related to the short-ranged nature of the base-stacking 
interactions |jj (see Sec. [IB). At even higher forces (/ > 80 pN) [J the DNA molecule becomes completely into the 
overstretched form with its folding angle peaked at (p = 0°. 

The force- induced axial fluctuations in DNA double- helix can be biologically significant. For example, it has been 
demonstrated that axial fluctuations in dsDNA enhance considerably the polymerization of RecA proteins along DNA 
chain [0J42M3]. An quantitative study on the coupling between RecA polymerization and DNA axial fluctuation is 



This threshold ft of over-stretch force is also consistent with a plain evaluation from base-stacking potential of e ~ ftro, 
ft ~ 90 pN. 







anticipated to be helpful. 

It seems that in the experiments fll],[l2| the transition to S-DNA occurs even more cooperatively and abruptly than 
predicted by the present theory (see Fig. ||). This may be related to the existence of single-stranded breaks (nicks) 
in the dsDNA molecules used in the experiments. Nicks in DNA backbones can lead to strand-separation or relative 
sliding of backbones [Tl|Jl^ ], and they can make the transition process more cooperative. However, the comprehensive 
agreement achieved in Figs. || and || indicates that such effects are only of limited significance. The elasticity of 
DNA is mainly determined by the competition between folding angle fluctuation and tangential fluctuation, which 
are governed, respectively, by the base-stacking interactions (e) and the axial bending rigidity («;*) in Eq. (jl3|). 

IV. ELASTIC PROPERTY OF SUPERCOILED DNA 

In the preceding section, we have discussed the elastic response of long DNA chains under the action of external 
forces. In the present section, we turn to study the elasticity of supercoiled DNA double-helix. For this purpose, in 
the experimental setup, all the possible nicks in the DNA nucleotide strands are ligated [Q, and a torque as well as 
an external pulling force is acting on one terminal of the DNA double-helix, which typically u ntwis ts or overtwists the 



original B-form double-helix to some extent and makes its total linking number (refer to Sec. II C for the definition of 
the linking number) less or greater than the equilibrium value. We say that such DNA molecules with deficit (excess) 
linking number are negatively (positively) supercoiled, and define the degree of supercoiling as 

Lk - Lk 

where Lko represents the linking number of a relaxed DNA of the same contour length. In living organisms, DNA 
molecules are often negatively supercoiled, with a linking number deficit of about a = —0.06. Thus, a detailed 
investigation on the mechanical property of supercoiled DNA molecules is not only of academic interest but can also 
help us to understand the possible biological advantages of negative supercoiling. 

A. Relationship between extension and supercoiling degree 

We focus on the property of highly extended DNA molecules whose tangent vectors fluctuate only slightly around 



the force direction Zo. According to what we have mentioned in Sec. [IC, in this case the approximate energy 
expression Eq. ( jl4| ) can be used. The external stretching force is restricted to be greater than 0.3 pN to make sure 
that the end-to-end distance of DNA chain approaches its contour length (see Fig. ||). Based on Eqs. ( fbi) ) and JA7| ), 
the Green equation for highly stretched and supercoiled dsDNA is then obtained to be 



dV(t,<p; s) 
ds 



d 2 d 2 fcostp p{ip) l p . 4 

t • z - — — - — £ sin <f 



U* p dt 2 At p dtp 2 k B T u k B T R 2 

r to r 2 . 2 

'M^f smip ~ 4k B Tt; ~d4> + iu* p {k B Tf sm ' 



*(t,W«). (23) 



where (6, 0) are the two directional angles of t as mentioned before in Sec. [II. Similar to what we have done in 



Sec. Ill , we can now express the above Green equation in matrix form using the combinations of spherical harmonics 
Yim{Q,<P) and as the base functions. The ground-state eigenvalue and eigenfunction of Eq. ( |23| ) can then be 

obtained numerically for given applied force and torque and the average extension be calculated through Eq. ( |l6| ) or 
through the following formula: 

(Z) = L / xo(t, <p)t ■ z cosy$ (t, <p)dtd<p, (24) 



where Xo(t,</?) is the g round-state left-eigenfunction of Eq. (p3|). ^ The writhing number Eq. ( jlj ) is calculated 



according to Eq. (A16) to be 



4 As remarked in Appendix |X| because the operator in the square brackets of Eq. (|2^) acting on ^(t, tp; s) is not Hermitian, the 
resulting matrix form of the operator may not be diagonalized by unitary matrices. Consequently, in general Xo(t, tp) 7^ $o(t) V 3 )- 



10 



(Wr) = L 16jrtkBT J xo(t, <p) ^ ^o(t, p)dtdp, (25) 

and average linking number is then calculated to be 

(Lk) = (Tw) + (Wr) = —J Xo sin <p$ dtd<p + ^~^f J Xo sin 2 9^ dtdcp. (26) 

Thus, after we have obtained the ground-state eigenvalue as well as its left- and right-eigenfunction numerically, we 
can calculate numerically all the quantities of our interest, for example, the average extension, the average supercoiling 
degree, the folding angle distribution (see also Appendix ^). The relation between extension and supercoiling degree 
can also be obtained by fixing the external force and changing the value of the applied torque. To calculate the 
ground-state eigenvalue and eigenfunctions of an asymmetric matrix turns out to be complicated and time-consuming. 
Fortunately, as we have calculated in Appendix each eigenfunction of Eq. (^3|) shares the same quantity to; the 
matrix for to = is still Hermitian and can be diagonalized by unitary matrices. The ground-state eigenvalue for 
to = is lower in several order than those for m ^ in the whole relevant region of external torque Y from — 5.0fcsT to 
5.0fcsT. Thus, actually we only need to consider the case of m = and in this case we still have xo(t, f) — Q^it, ip). 



The whole procedure we have performed in Sec. Ill can safely be repeated in this section, and the relationships 



between force and extension and between torque and linking number can be consequently calculated. 
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FIG. 5. Extension vs supercoiling relations at fixed pulling forces for torsionally constrained DNA. The parameters for the 
curves are the same as Fig. [| and experimental data is from Fig. 3 of |l3| (symbols) . 

To make the calculation further easier and also to make sure that the above-mentioned calculation is indeed correct, 
we introduce here an approximate method which reduces the computational complexity considerably. It turns out 
that the calculated results using this method are in considerable agreement with the above-mentioned precise method. 
In the experiment of Rcf. [OL the applied external forces change in the region of 0.3 pN to 10 pN. In this region, as 
demonstrated in Figs. || andf|, both the tangential (t) and the folding angle (up) fluctuations of dsDNA molecules are 
small. Taking into account this fact, then in Eq. (|14|), the energy term K*(dt/ds) 2 can be approximately calculated 
to be K*((dt x /ds) 2 + (dt v /ds) 2 ), and ft ■ z cosip ~ fcosip — /(cosy)(t| + t 2 y )/2. Thus, Eq. ( |l4j) is decomposed 
into two "independent" parts. The first part is only related to (p. At each value of / and T, we can calculate the 
average quantities (cost/?) and (sinip) based on this energy using the method of path integral. The second part is 
quadratic in t x and t y , therefore the average values of t • z and Wr (Eq. (|ll])) can be obtained analytically. Using 
this decomposition and preaveraging technique, the average extension and average supercoiling degree can both be 
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calculated at each value of external force and torque, and the relation between extension and linking number at fixed 
forces can be then obtained. 

The theoretical relationship between extension and supercoiling degree is shown in Fig. and compared with the 
experiment of Strick et al. [113] . In obtaining these curves, the values of the parameters are the same as those used 
in Fig. U and no adjustment has been done to fit the experimental data. We find that in the case of negatively 
supercoiled DNA, the theoretical and experimental results are in quantitative agreement, indicating that the present 
model is capable of explaining the elasticity of negatively supercoiled DNA; in the case of positively supercoiled DNA, 
the agreement between theory and experiment is not so good, especially when the external force is relatively large. In 
our present work, we have not considered the possible deformations of the nucleotide basepairs. While this assumption 
might be reasonable in the negatively supercoiled case, it may fail for positively supercoiled DNA chain, especially at 
large stretching forces. The work done by Allemand et al. [15) suggested that positive supercoiled and highly extended 
DNA molecule can take on Pauling-like configurations with exposed bases. To better understand the elastic property 
of positively supercoiled DNA, it is certainly necessary for us to take into account the deformations of basepairs. 

For negatively supercoiled DNA molecule, both theory and experiment reveal the following elastic aspects: (a) 
When external force is small, DNA molecule can shake off its torsional stress by writhing its central axis, which can 
lead to an increase in the negative writhing number and hence restore the local folding manner of DNA strands to 
that of B-form DNA: (b) However, writhing of the central axis causes shortening of DNA end-to-end extension, which 
becomes more and more unfavorable as the external force is increased. Therefore, at large forces, the torsional stress 
caused by negative torque (supercoiling degree) begins to unwind the B-form double-helix and triggers the transition 
of DNA internal structure, where a continuously increasing portion of DNA takes on some certain new configuration 
as supercoiling increases, while its total extension keeps almost invariant. Our Monte Carlo simulations have also 
confirmed the above insight E3 . 




FIG. 6. Folding angle distributions for negatively supercoiled DNA molecule pulled with a force of 1.3 pN. 

What is the new configuration? According to the opinion of Ref. [ fL3f , such new configuration corresponds to 
denatured DNA segments, i.e., negative torque leads to breakage of hydrogen bonds between the complementary 
DNA bases and consequently to strand-separation. They have also done an elegant experiment in which short single- 
stranded homologous DNA segments are inserted into the experimental buffer frL4[] . They found that, in confirmative 
with their insight, these homologous DNA probes indeed bind onto negatively supercoiled dsDNA molecules. Recently, 
Leger et al. [lrj| ] have also done experiment on single dsDNA molecules. To explain qualitatively their experimental 
result, they found that left-handed Z-form DNA should be considered as a possible configuration for negatively 
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supercoiled dsDNA chain, while the molecules need not be denatured. As seen in Fig. H, although the present 
model has not taken into account the possibility of strand-separation, it can quantitatively explain the behavior of 
negatively supercoiled DNA. Therefore, it may be helpful for us to investigate the possibility of formation of left- 
handed configurations based on our present model. This effort is done in the next subsection, where the energetics of 
such configurations will also be discussed. 



B. Possible left-handed DNA configurations 



We have mentioned in Sec. [I A that left-handed configurations correspond to if < in our present mo del (s ee also 
Fig. [I]). Based on the present model, then information about the new configuration mentioned in Sec. fV A| can be 
revealed by the folding angle distributions P{(p), 
calculated as follows: 



as have been discussed in Sec. [II. In the present case, P(y>) is 



P(<P) = / Xo(t,fp)*o(t,v)dt, 



(27) 



where, as mentioned in Sec. IV A, 3>o(t, an d Xo(t, f) are, respectively, the ground-state right- and left-eigenfunction 
of Eq. (HH); and actually, xo(t, <p) = $o(t, <p). 

The calculated folding angle distribution El is shown in Fig. ||, which is radically different with that of torsionally 
relaxed dsDNA molecules shown in Fig. ||. It has the following aspects: When the torsional stress is small (with the 
supercoiling degree \a\ < 0.025), the distribution has only one narrow and steep peak at <p ~ +57.0°, indicating that 
DNA is completely in B-form. With the increase of torsional stress, however, another peak appears at <p ~ —48.6° 
and the total probability for the folding angle to be negatively-valued increases gradually with supercoiling. Since 
negative folding angles correspond to left-handed configurations, the present model suggests that, with the increasing 
of supercoiling, left-handed DNA conformation is nucleated and it then elongates along the DNA chain as B-DNA 
disappears gradually. The whole chain becomes completely left-handed at a ~ —1.85. 




External torque (unit: k B T) External torque (unit: k B T) 



FIG. 7. (A), The sum of the average base-stacking and torsional energy per basepair at force 1.3 pN. For highly extended 
DNA only these two interactions are sensitive with torque. (B), The relation between DNA supercoiling degree and external 
torque at force 1.3 pN. 

It is worth to be noticed that, (a) as the supercoiling degree changes, the positions of the two peaks of the folding 
angle distribution remain almost fixed and, (b) between these two peaks, there exists an extended region of folding 
angle from to 7r/6 which always has only extremely small probability of occurrence. Thus, a negatively supercoiled 
DNA can have two possible stable configurations, a right-handed B-form and a left-handed configuration with an 
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average folding angle ~ —48.6°. A transition between these two structures for a DNA segment will generally lead to 
an abrupt and finite variation in the folding angle. 

To obtain the energetics of such transitions, we have calculated how the sum of the base-stacking energy and 
torsional energy, (k/R 2 ) sin 4 tp + p{ip) — (T/R) siny, changes with external torque fig] . Figure ^a shows the numerical 
result, and Fig. [7p demonstrates the relation between supercoiling degree and external torque. (In both figures the 
external force is fixed at 1.3 pN.) From these figures we can infer that, (a) for negative torque less than the critical 
value r c ~ —3.8 k^T, DNA can only stay in B-form state; (b) near this critical torque, DNA can either be right- or 
be left-handed and, as negative supercoiling increases (see Fig. 0b) more and more DNA segments will stay in the 
left-handed form, which is much lower in energy (~ —2.0 k^T per basepair) but stable only when torque reaches F c ; 
(c) for negative torque greater than T c DNA is completely left-handed. 

Nevertheless, we should emphasize that the above calculations are all based on our present model which has assumed 
that nucleotide basepairs do not break. Figure ^ indicates that for negatively supercoiled DNA chains, the extension 
vs supercoiling degree relation can be quantitatively explained by the present model; and Fig. ^ reveals the reason 
of the quantitative agreement is that the present model allows the possibility of occurrence of left-handed DNA 
configurations. At the present time, to say that negatively supercoiled DNA will prefer left-handed configurations 
rather than denaturation and strand-separation is premature. To clarify this question, the present model should be 
improved to consider the deformations of DNA basepairs. In the experimental side, it might also be helpful to measure 
precisely the critical torque at which the elastic behavior of negatively supercoiled DNA changes abruptly and compare 
the measured results with the value calculated in the present work. (In the earlier experiment of Allemand et al. |fL5| , 
the critical torque is estimated to be ~ —2ksT by assuming that the torsional rigidity of dsDNA to be 75 nm and 
that torsional stress builds up linearly along the DNA chain. This value, however, maybe not precise enough since the 
torsional rigidity of dsDNA molecule is not a precisely determined quantity and the values giving by different groups 
are scattered widely.^] 

The structural parameters of the left-handed configuration suggested by Figs. || and]?] are listed in table | Ej| and 
compared with those of Z-form DNA Q . The strong similarity in these parameters suggests that the torque-induced 
left-handed configurations, if they really exist, belong to Z-form DNA p ]l6| Jt5[l . 



TABLE I. For the torque-induced left-handed DNA configuration, the average rise per basepair (d), the pitch per turn of 
helix, and the number of basepairs per turn of helix {num) are calculated and listed under different external forces and torques. 
The last row contains the corresponding values for Z-form DNA. 



force (pN) 


torque (fcsT) 


d(A) 


pitch (A) 


num 


1.3 


-5.0 


3.59 


41.20 


11.48 


1.0 


-5.0 


3.57 


40.93 


11.44 


1.3 


-4.0 


3.83 


46.76 


12.19 


1.0 


-4.0 


3.82 


46.38 


12.15 


Z-form: 




3.8 


45.6 


12 



5 Another possibility may be that we have over-estimated the value of the critical torque. It maybe possible that the transition 
from right- to left-handed configurations is initiated in the weaker AT-rich regions, whose value of e should be less than the 
average value taken by the present paper. 
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V. CONCLUSION 



In this article, we have presented an elastic model for double-stranded biopolymers such as DNA molecules. The 
key progress is that the bending deformations of the backbones of DNA molecules and the base-stacking interactions 
existing between adjacent DNA basepairs are quantitatively considered in this model, with the introduction of a 
new structural parameter, the folding angle ip. This model has also qualitatively taken into account the effects of 
the steric effects of DNA basepairs and electrostatic interactions along DNA chain. In calculation technique, the 
model is investigated using path integral method; and Green equations similar in form to the Schrodinger equation in 
quantum mechanics are derived and their ground-state eigenvalues and eigenfunctions obtained by precise numerical 
calculations. The force-extension relationship in torsionally relaxed and the extension-linking number relationship in 
torsionally constrained DNA chains are studied and compared with experimental results. This work demonstrated that 
DNA molecule's entropic elasticity and highly extensibility, as well as the elastic property of negatively supercoiled 
DNA can all be quantitatively explained by the present theory. The comprehensive agreement between theory and 
experiments indicated that the short-ranged base-stacking interactions are very important in determining the elastic 
response of double-stranded DNA molecules. The present work showed that highly extended and negatively supercoiled 
DNA molecules can be left-handed, probably in the Z-form configuration. A possible way to check the validity of 
this opinion is to measure the critical external torque at which the transition between B-form DNA and the new 
configuration takes place. 

The present work regarded DNA basepairs as rigid objects and did not consider their possible deformations and the 
possibility of strand-separation. The comprehensive agreement between theory and experiments indicates that this 



approach is well justified in many cases. However, as already mentioned in Sec. II A. under some extreme conditions, 
this assumption may not be appropriate. For example, recent experimental work of Allemand et al. [ [l5| showed that 
positively supercoiled DNA under high applied force can take on Pauling-like configurations with exposed bases. In 
this case the basepairing of DNA is severely distorted, and because the present model has not taken into account 
the possible deformations of the basepairs, the theoretical results on positively supercoiled DNA molecules were not 
in quantitative agreement with experiment (see Fig. ^J). Furthermore, although the present work showed that left- 
handed DNA configurations can be stabilized by negative torques, much theoretical work is still needed to calculate 
the denaturation free energy and be compared with the free energy of left-handed DNA configurations. 
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APPENDIX A: PATH INTEGRAL METHOD IN POLYMER PHYSICS 

In this appendix we review some basic ideas on the application of path integral method to the study of polymeric 
systems p^] . Consider a polymeric string, and suppose its total "arclength" is L, and along each arclength point s 
one can define a n-dimcnsional "vector" r(s) to describe the polymer's local state at this point.0 We further assume 
that the energy density (per unit arclength) of the polymer can be written as the following general form: 

p e (r,.s) = |(J) 2 + A(r). f g +V(r), (Al) 

where T^(r) is a scalar field and A(r) is a vectorial field. The total partition function of the system is expressed by 
the following integration: 



B(L)= / dr / ^ / (r / )G(r / ,L;r i ,0)&(r i )dr i , (A2) 



6 For example, in the case of a flexible Gaussian chain, r is a three-dimensional position vector; in the case of a semiflexible 
chain such as the wormlike chain , r is the unit tangent vector of the polymer and is two-dimensional. 
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where 4>i(r) and 4>f{r) are, respectively, the probability distributions of the vector r at the initial (s = 0) and final 
(s — L) arclength point; G(r, s; r', s') is called the Green function, it is defined in the following way: 



G(r, s; r', s') = £ 2>[r»] exp[-p £ ds" p e (v" , S ")], 



(A3) 



where integration is carried over all possible configurations of r", and j3 = l/fc^T is the Boltzmann coefficient. It can 
be verified that the Green function defined above satisfies the following relation |39] : 



G(r, s; r', s') = J dr"G(r, s; r", s")G(r", s"; r', s'), (s' < s" < s). 



The total free energy of the system is then expressed as 

T = —ksT In S. 



(A4) 



(A5) 



To calculate the total partition function S, we define an auxiliary function ^(r, s) and call it the wave function 
because of its similarity with the true wave function of quantum systems. Suppose the value of ^ at arclength point 
s is related to its value at s' through the following formula such thatQ 



*(r, s) = J dr'G(r, s; r', s')*(r', s'), (s > s') 



then, we can derive from Eqs. (Al), (|A3|), and (A6) that B9]: 



d~s 



2m/3 



pV(v) 



A(r)-V r V r -A(r) /3A 2 (r) 



2m 



2m 



\l/(r, s) = H^>(r, s) 



(A6) 



(A7) 



Equation (A7) is called the Green equation, it is very similar with the Schrodinger equation of quantum mechanics 



MO). However, there is an important difference. In the case of A(r) ^ 0, the operator H in Eq. (A7) is not Hcrmitian. 
Therefore, in this case the matrix form of the operator H may no t be diagonalized by unitary matrix. 

Denote the eigenvalues and the right-eigenfunctions of Eq. ( |A7| ) as — gt and \i) = $i(r) (i = 0, 1, ■ ■ •), respectively. 
Then it is easy to know, from the approach of quantum mechanics, that 



-gi(s-s') 



(A8) 



where (i\ = Xi( r ) — 0, 1, • • •) denote the left-eigenfunctions of Eq. (A7), which satisfy the following relation: 

(i\i')= [ drxi(r)®i'(r)=5i'. 



In the case where H is Hcrmitian (i.e., A(r) = 0), then we can conclude that 

Xi(r) =$?(r). 



From Eqs. (|A2|), dA6|), (|A7|), and (JA8|) we know that 

S(L) = J dv J dT'G{T,L ] v\s')ct> { {v)cj> i {v')=Y J ^f\i)^i)e 

= e-^^flO)^) (for L » l/(g x -g )). 
Consequently, for long polymer chains the total free energy density is just expressed as 

TjL = k B Tg , 



-9iL 



(A9) 
(A10) 



7 In fact, the choice of the wave function ^(r, s) is not limited. Any function determined by an integration of the form of Eq. 
(Kq) can be viewed as a wave function. 
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and any quantity of interest can then be calculated by differentiation of T . For example, the average extension of a 
polymer under external force field / can be calculated as (Z) = dF/df = LksTdgo/ df . 

We continue to discuss another very important quantity, the distribution probability of r at arclength s, P(r,s). 
This probability is calculated from the following expression: 



P(r,s) = 



J dv f J drj(j)f(rf)G(rf,L;r,s)G(r,s;ri,0)(j)i(ri) 
I dv f I d n4>f(vf)G(r f ,L; r u 0)^0,) 



(All) 



Based on Eqs. (A6) and (A8) we can rewrite Eq. ( All ) in the following form: 

Jdr f J drt J dr'c/> f (r f )G(r f ,L; r', s)S(r' - r)G(r, s; r i; 0)^ 



P(r,s) = 



I dr f I dr i0/( r /)G(r/, L; r u 0)^(r 4 ) 
EEf^/hX^l^lXmW^W exp[- 3m (i - s) -g n s] 



Y,{<t>f \ m ) ( m \4>i) exp(-g m L) 



(A12) 



For the most important case of < s < L, Eq. (|A12| ) then gives that the probability distribution of r is independent 
of arclength s, i.e., 



P{r, s) = xo(r)$o(r) (for < s < L). 



(A13) 



With the help of Eq. ( A13), the average value of a quantity which is a function of r can be obtained. For example, 



(Q(s)) = / drQ(r)P(r,s) = / dr X o(r)Q(r)<I>o(r) = <0|Q|0) 



and 



Q(r(s))ds) 



(Q(s))ds = L(0\Q\0) (for L > l/( 9l - go)). 



(A14) 



(A15) 



Finally, we list the formula for calculating (B(r) • dr/ds), here B(r) is a given vectorial field. The formula reads: 

1 



r ■ (HB) - B • (Hr) |0) (for L » l/( 9l - go)). 



(A16) 



APPENDIX B: MATRIX FORMALISM OF THE GREEN EQUATION EQ. (23) 



In Eq. ( |23| ) (and also Eq. (|15|)) the variables t and ip are coupled together. Denote the operator in the square 
brackets of Eq. ( |23] ) as H. We express this operator in matrix form. To this end we choose the base functions of this 
system to be the combinations of spherical harmonics Y[ m (9, 4>) and f n (<p) (see Eq. (po|)). 

For m = 0, the base functions are expressed to be 

|t> = |2 ■ N v + n> = |/; n) = Y l0 (6, <j>)U(<p)\ (Bl) 

and for m = 1, 2, • ■ •, the base functions are expressed to be 

|t> = \2(l - m) ■ N v + 2(n - 1) + k) = \l; n; k) = y£\o, <f>)f n ((f>). (B2) 

In the above two equations, n — 1,2, ■ ■ ■ , N v with N v set to be 60 in our present calculations; / = 0, 1, • • • , Ni — 1 for 
the case of m = (with Ni set to 30), or I = m, m + 1, • • • , m + JVj — 1 for the case of m / (with Ni set to 15); 
k = 1,2; and 

1 Pi (cos ( 




21 + 1 (l-m)\ Dmf . , 

— C, if (cos 9) sm m</> . 

2ir [I + my. 
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With the above definitions, we then obtain that, for to = 0: 



(i P \(-H)\i) = (l p ;n p \(-H)\l;n) 



f 
k B T 

S K n P 



1(1 + 1) n 2 TT 2 



4/'f 



U p a 2 



S l , +1 aio + 5 l , 1 a;_i. 

lr> ' IT) 



(n p \ sin <p\n) 



*-p 4 - 
^cos <p 



r 



82 



"Rk B T 

r 2 



k B T 
n p \ cos (p\n) 

<t(l 



and for m =/= 0: 

(ip|(-F)|i) 



» 16£;(fc B T) 



(l p ;n p ;k p \(-H)\l;n;k) 



H,0 



a l-l o 



(B3) 



f 



1(1 + 1) 

41; 



A£ p a 2 



k B T 
+ S l A(np\ 

r 



p(<p) I 



( rip) shnp\n) 



k B T 



P 4 ~ 
-cos <p 



\n) 



5" 5j£ 



Rk B T 
T 2 



(n p \ cos (p\n) + 5\ 51 



' 4^;fc B T 



(k — k p )m 



" 16£;(k B T) 2 



*L(l-a?, 



i-2. 



(B4) 



In Eqs. (B3) and (B4), the following notations are used: 



aur, 



(n p \y((f)\n 



(l + l) 2 -m 2 
(2/ + l)(2Z + 3)' 

If 

a Jo 



sm(- )y(<f) sm( ), 



where y((p) is any function of (p. 

The ground-state eigenvalues of the matrices Eq. (B3) and Eq. (B4) have been calculated at force 1.3 pN in the 
whole relevant region of T from —5.0k B T to +5.0k B T. We have found that the ground-state eigenvalues for the case 
of m = are of the order of —10 nm" 1 ; while those for the case of to ^ are of the order of 10 6 nm _1 (data not 
shown). Thus, we can conclude with confidence that the "low energy" eigenstates of the system all have the same 
to = 0. This can greatly reduce the calculation tasks. 

As a final point, here we demonstrate how to calculate the average values of the quantities of our interest. Suppose 
y(t,(f) is a quantity whose average value we are interested in. Denote U as the unitary matrix which can diagonalizc 
the matrix Eq. (B3). (Because the matrix Eq. (B3) is real symmetric, U is actually a real orthogonal matrix.) Then 
its column vector u(i) — J2j i) \j) (i = 1, 2, • • •) corresponds to the i-ih eigenvector of Eq. (|B3]). Consequently, 



we can calculate based on Eq. (A14) that 

(y(t, $)) =Y, U( ?p> iXvIK*. <f)\i)- 



(B5) 
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